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The QCD equation of state at finite T and /i 
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We calculate the pressure (p), the energy density (e) and the baryon density (ns) of QCD at finite temperatures 
(T) and chemical potentials (fi). The recently proposed overlap improving multi-parameter reweighting technique 
is used to determine observables at nonvanishing chemical potentials. Our results are obtained by studying 
Uf =2+1 dynamical staggered quarks with semi-realistic masses on A'^t = 4 lattices. 
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Figure 1. The best weight lines on the /i-/3 plane. 
In the middle we indicate the transition line. Its 
first dotted part is the crossover region. The blob 
represents the critical endpoint, after which the 
transition is of first order. The integration paths 
used to calculate p are shown by the arrows along 
the f3 axis and the best weight lines. 

QCD at nonzero density is easily formulated on 
the lattice by multiplying the forward/backward 
links by exp(±/^). However, standard importance 
sampling based Monte-Carlo techniques can not 
be used at fi y^O. Up to now, no technique was 
suggested capable of giving the equation of state 
(EOS) at /i^O, which is essential to describe the 
quark gluon plasma (QGP) formation at heavy 
ion collider experiments. Results are only avail- 
able for ^=0 (e.g. |l|,||]) at T^^O. 

The overlap improving multi-parameter 
reweighting B opened the possibility to study 




Figure 2. p normalised by as a function of 
T/Tc at /i = The iVf=4 SB hmit is also shown. 

lattice QCD at nonzero T and fi. First one pro- 
duces an ensemble of QCD configurations at /i=0 
and T^O. Then reweighting factors are deter- 
mined at /i 7^ and at a lowered T. The idea is 
expressed in terms of the partition function 



Z(/3,/i, m) 



J VU exp[-Sg{P)]det M{fi,m) 
PC/exp[-5g(/3o)]detA/(Ai = 0,m) (1) 



eM-Sg{P) + SgiPo)]- 



det Af(/i, to) 
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where to is the quark mass, Sg is the action of 
the gluonic field (C/), at /? = coupling. At 
nonzero one gets a complex fermion determi- 
nant det M which spoils importance sampling. 
Thus, the first line of eq. (|^), /i t^O, is rewritten 
in a way that the second line of eq. Im) is used as 
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Figure 3. Ap = ^ 0, T) - = 0, T) nor- 
malised by T"^ as a function of T/Tc for /is^lOO, 
210, 330, 410 MeV and /ii3=530 MeV (from bot- 
tom to top). 

an integration measure (at ii=0 with importance 
sampling) and the curly bracket is measured on 
each independent configuration and is interpreted 
as a weight factor {w{P, ^,m,U)}. In order to 
maximise the accuracy of Z the reweighting is 
performed along the best weight lines on the /z- 
(3 plane. These best weight lines are defined by 
minimising the spread of log w. 

Using the above technique, transition (or 
hadronic/QGP) configurations are reweighted to 
transition (or hadronic/QGP) configurations as 
illustrated by Fig. ^ The technique works for 
temperatures at, below and above the transition 
temperature (Tc). By using the reweighting tech- 
nique, the phase diagram Q and the location of 
the critical endpoint j|] was given. Using a Taylor 
expansion around /i=0, T^^O for small fj, can be 
used to determine thermal properties Q. A dif- 
ferent method, analytic continuation from imagi- 
nary fi, confirmed results of |^,|| on the /x-T dia- 
gram 0. 

We use 4 -N^ lattices at T^O with iV^=8,10,12 
for reweighting and we extrapolate to V^oo us- 
ing the available volumes (V). At T=0 lattices 
of 24 • 14"^ are taken for vacuum subtraction and 
to connect lattice parameters to physical quanti- 
ties. 14 different f3 values are used, which corre- 
spond to T/Tc = 0.8, . . . , 3. Our T=0 simulations 
provided Rq and a. The lattice spacing at Tc is 
«0. 25-0. 30 fm. We use 2+1 flavours of dynami- 
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Figure 4. Ap of the QCD plasma normalised by 
Ap of the free gas (SB) as a function of T/Tc for 
the same /is values as in Fig. 3. 

cal staggered quarks. While varying (3 (thus the 
temperature) we keep the physical quark masses 
constant at niud ~ 65 MeV and iris ~ 135 MeV 
(the pion to rho mass ratio is mj^/rrip RiO.66). 

The determination of the equation of state 
at fj, 7^0 needs several observables O at non- 
vanishing fi values. This can be calculated by 
using the weights of eq. (|l|) 



0(/3,Ai,m) 



J:{w{p,^i,m,U)}Oif3,^l,m,U) 
J2{w{P,fJ.,m,U)} 



(2) 



p can be obtained from the partition func- 
tion as p=T ■ d\ogZ/dV which can be written 
as p—{T/V) ■ \ogZ for large homogeneous sys- 
tems. On the lattice we can only determine 
the derivatives of log^ with respect to the pa- 
rameters of the action (/3, m, /i). Using the fol- 
lowing notation (0(/3, /i, m))= 0(/3, /x, TO)j,_^g — 
0(/3, /X = 0, rnjrp^f^. p can be written as an inte- 
gral §: 



P_ 



1 



TW 
d{\ogZ) 
dp 



d{f3,m,fi) 

djlogZ) 
dm 



(3) 



9(log^) 
9/1 



The integral is by definition independent of the 
integration path. The chosen integration paths 
are shown on Fig ^ 

The energy density can be written as e = 
(T^/V) ■ d{log Z) jdT + {[iTlY] ■ d{\og Z) Idp,. By 
changing the lattice spacing T and V are simul- 
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Figure 5. (e - 3p)/T4 at ^b=0, 210, 410 MeV 
and 530 MeV versus T/Tc (from bottom to top). 

taneously varied. The special combination e — 3p 
contains only derivatives with respect to a and fi: 



d\og{Z) 



TW da 



fi a(logZ) 



•(4) 



The quark number density is n = (T/V) ■ 
d\og{Z) / dfji which can be measured directly or 
obtained from p (baryon density is ns—n/S and 
baryonic chemical potential is /iB=3/i). 

We present direct lattice results on p{^ — 0,T), 
Ap(/x,r) = p{p ^ 0,T)-pif, = 0,r), e{fi,T)- 
3p{iJ,,T) and nB(^,r). Note, that in [§ addi- 
tional overall factors were used to help the phe- 
nomenological interpretation. Our statistical er- 
rorbars are rather small, sometimes smaller than 
the thickness of the lines. 

Fig. ^ shows p at /i=0. On Fig. || we present 
Ap/T'^ for five different /Lt values. Fig. ^ gives 
Ap{n,T/Tc) normalised by Ap^^ = Ap{^,T ^ 
oo). Notice the interesting scaling behaviour. 
Ap/ Ap^^ depends only on T and it is practically 
independent of /j, in the analysed region. Fig. |^ 
shows e-3p normalised by T'*, which tends to zero 
for large T. Fig. || gives the baryonic density as 
a function of T/Tc for different /i-s. As it can be 
seen the densities exceed the nuclear density by 
up to an order of magnitude. 

An important finding concerns the applicability 
of our reweighting method: the maximal /i scales 
with the volume as /Xmax • a ~ {Nt ■ TV^ )~°-^^. If 
this behaviour persists, one could -in principle- 
approach the true continuum limit (a ^ l/iV* 
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Figure 6. ub/T^ versus T/Tc for the same 
values as in Fig. 3 (from bottom to top). 

{Nt ■ iV3)-0-25, thus Aimax«COnSt.). 

Future analyses should be done at smaller lat- 
tice spacings and quark masses. A detailed ver- 
sion of this work can be found elsewhere ^ . 
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